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ABSTRACT 


We use 32 age measurements of passively evolving galaxies as a function of 
redshift to test and compare the standard model (ACDM) with the Rh = ct 
Universe. We show that the latter fits the data with a reduced y^ of = 0.435 for 
a Hubble constant H 0 = 67.21 4 q km s -1 Mpc -1 . By comparison, the optimal 
flat ACDM model, with two free parameters (including Q m = 0.121 q n and Hq = 
94.3^3^ km s -1 Mpc -1 ), fits the ag e-z data with a reduced x\ Q i — 0.428. Based 
solely on their y 4of values, both models appear to account for the data very 
well, though the optimized ACDM parameters are only marginally consistent 
with those of the concordance model (D m = 0.27 and H 0 = 70 km s -1 Mpc -1 ). 
Fitting the age -2 data with the latter results in a reduced xI of = 0.523. However, 
because of the different number of free parameters in these models, selection tools, 
such as the Akaike, Kullback and Bayes Information Criteria, favour R.y = ct over 
ACDM with a likelihood of ~ 66.5% — 80.5% versus ~ 19.5% — 33.5%. These 
results are suggestive, though not yet compelling, given the current limited galaxy 
age -2 sample. We carry out Monte Carlo simulations based on these current age 
measurements to estimate how large the sample would have to be in order to 
rule out either model at a ~ 99.7% confidence level. We find that if the real 
cosmology is ACDM, a sample of ~ 45 galaxy ages would be sufficient to rule 
out Ry, = ct at this level of accuracy, while ~ 350 galaxy ages would be required 
to rule out ACDM if the real Universe were instead Ry, = ct. This difference in 
required sample size reflects the greater number of free parameters available to 
fit the data with ACDM. 


Subject headings: cosmology: theory, observations - early universe - galaxy: general 
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1. Introduction 

In recent years, the cosmic evolution has been studied using a diversity of observational 
data, including cosmic chronometers (Melia & Maier 2013), Gamma-ray bursts (GRBs; 
Wei et al. 2013), high -z quasars (Melia 2013, 2014a), strong gravitational lenses (Wei 
et al. 2014; Melia et al. 2015), and Type la SNe (Wei et al. 2015). In particular, the 
predictions of ACDM have been compared with those of a cosmology we refer to as the 
R h = ct Universe (Melia 2007; Melia & Shevchuk 2012). In all such one-on-one comparisons 
completed thus far, model selection tools show that the data favour R h = ct over ACDM. 

The R h = ct Universe is a Friedmann-Robertson-Walker (FRW) cosmology that has 
much in common with ACDM, but includes an additional ingredient motivated by several 
theoretical and observational arguments (Melia 2007; Melia & Abdelqader 2009; Melia & 
Shevchuk 2012; see also Melia 2012a for a more pedagogical treatment). Like ACDM, it 
adopts an equation of state p = wp, with p = p m + p r + p de (for matter, radiation, and 
dark energy, respectively) and p = p m + p r + p de , but goes one step further by specifying 
that w = (p r /3 + WdePde)/P = —1/3 at all times. Here, p is the total pressure and p is the 
total energy density. One might come away with the impression that this equation of state 
cannot be consistent with that (i.e., w = [p r /3 — Pa]/p) in the standard model. But in fact 
if we ignore the constraint w = —1/3 and instead proceed to optimize the parameters in 
ACDM by fitting the data, the resultant value of w averaged over a Hubble time is actually 
—1/3 within the measurement errors (Melia 2007; Melia & Shevchuk 2012). In other words, 
though w = (p r /3 — Pa)/p hi ACDM may be different from —1/3 from one moment to the 
next, its value averaged over the age of the Universe equals what it would have been in 
i?h = ct all along (Melia 2015). 

In this paper, we continue to compare the predictions of ACDM with those in the 
i?h = ct LIniverse, this time focusing on the age-redshift relationship, which differs from 



one expansion scenario to another. Though the current age of the universe may be similar 
in these two cosmologies, the age versus redshift relationship is not, particularly at high 
redshifts (Melia 2013, 2014b). Previous work with the age estimates of distant objects has 
already provided effective constraints on cosmological parameters (see, e.g., Alcaniz & Lima 
1999; Lima & Alcaniz 2000; Jimenez & Loeb 2002; Jimenez et al. 2003; Capozzicllo et al. 
2004; Friaga et ah 2005; Simon et ah 2005; Jain & Dev 2006; Pires et ah 2006; Dantas 
et ah 2007, 2009, 2011; Samushia et ah 2010). For example, using the simple criterion 
that the age of the Universe at any given redshift should always be greater than or equal to 
the age of the oldest object(s) at that redshift, the measured ages were used to constrain 
parameters in the standard model (Alcaniz & Lima 1999; Lima & Alcaniz 2000; Jain & Dev 
2006; Pires et ah 2006; Dantas et ah 2007, 2011). Cosmological parameters have also been 
constrained by the measurement of the differential age Az/At, where Az is the redshift 
separation between two passively evolving galaxies having an age difference At (Jimenez 
& Loeb 2002; Jimenez et ah 2003). And the lookback time versus redshift measurements 
for galaxy clusters and passively evolving galaxies have been used to constrain dark energy 
models (Capozzicllo et ah 2004; Simon et ah 2005; Dantas et ah 2009; Samushia et 
ah 2010). This kind of analysis is therefore particularly interesting and complementary 
to those mentioned earlier, which are essentially based on distance measurements to a 
particular class of objects or physical rulers (see Jimenez & Loeb 2002, for a discussion on 
cosmological tests based on relative galaxy ages). 

Age measurements of high -2 objects have been valuable in constraining the cosmological 
parameters of the standard model even before dark energy was recognized as an essential 
component of the cosmic fluid (see, e.g., Bolte & Hogan 1995; Krauss & Turner 1995; 
Dunlop et ah 1996; Alcaniz & Lima 1999; Jimenez & Loeb 2002). Of direct relevance to the 
principal aim of this paper is the fact that, although the distance-redshift relationship is 
very similar in ACDM and the A h = ct Universe (even out to z > 6 — 7; Melia 2012b; Wei et 
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al. 2013), the age-redshift dependence is not. This difference is especially noticeable in how 
we interpret the formation of structure in the early Universe. For example, the emergence 
of quasars at z > 6, which are now known to be accreting at, or near, their Eddington limit 
(see, e.g., Willott et al. 2010; De Rosa et al. 2011). This presents a problem for ACDM 
because it is difficult to understand how ~ 10 9 M 0 supermassive black holes could have 
appeared only 700-900 Myr after the big bang. Instead, in R h = ct , their emergence at 
redshift ~ 6 corresponds to a cosmic age of > 1.6 Gyr, which was enough time for them 
to begin growing from ~ 5 — 20 M 0 seeds (presumably the remnants of Pop II and III 
supernovae) at z < 15 (i.e., after the onset of re-ionization) and still reach a billion solar 
masses by z ~ 6 via standard, Eddington-limited accretion (Mclia 2013). 

In this paper, we will broaden the base of support for this cosmic probe by 
demonstrating its usefulness in testing competing cosmological models. Following the 
methodology presented in Dantas et al. (2011), we will use 32 age measurements of 
passively evolving galaxies as a function of redshift (in the range 0.117 < z < 1.845) to test 
the predicted age-redshift relationship of each model. From an observational viewpoint, 
because the age of a galaxy must be younger than the age of the Universe at any given 
redshift, there must be an incubation time, or delay factor r, for the galaxy to form after 
the big bang. In principle, there could be a different T t for each object i since galaxies can 
form at different epochs. However, the simplest approach we can take is to begin with the 
assumption made in earlier work (see, e.g., Dantas et al. 2009, 2011; Samushia et al. 2010), 
i.e., we will adopt an average delay factor (r) and use it uniformally for every galaxy. But 
we shall also consider cases in which the delay factors r t are distributed, and study the 
impact of this non-uniformity on the overall fits to the data. 

We will demonstrate that the current sample of galaxy ages favours the R,\, = ct 
Universe with a likelihood of ~ 66.5 — 80.5% of being correct, versus ~ 19.5 — 33.5% 
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for ACDM. Though this result is still only marginal, it nonetheless calls for a significant 
increase in the sample of suitable galaxy ages in order to carry out more sophisticated 
and higher precision measurements. We will therefore also construct mock catalogs to 
investigate how big the sample of measured galaxy ages has to be in order to rule out one 
(or more) of these models. 

The outline of this paper is as follows. In § 2, we will briefly describe the age-redshift 
test, and then constrain the cosmological parameters—both in the context of ACDM and 
the R h = ct universe, first using a uniform (average) delay factor (r) (§ 3), and then a 
distribution of r* values (§ 4). In § 5, we will discuss the model selection tools we use to test 
ACDM and the Rh = ct cosmologies. In § 6, we will estimate the sample size required from 
future age measurements to reach likelihoods of ~ 99.7% and ~ 0.3% (i.e., 3a confidence 
limits) when using model selection tools to compare these two models, and we will end with 
our conclusions in 5 7. 


2. The age-redshift test 


The theoretical age of an object at redshift z is given as 

3 dzJ 


t (2,p) = 




( 1 ) 


where p stands for all the parameters of the cosmological model under consideration and 
H{z, p) is the Hubble parameter at redshift z. From an observational viewpoint, the total 
age of a given object (e.g., a galaxy) at redshift z is given by t ohs (z) = tc{z) + r, where 
tc(z) is the estimated age of its oldest stellar population and r is the incubation time, or 
delay factor, which accounts for our ignorance about the amount of time elapsed since the 
big bang to the epoch of star creation. 


To compute model predictions for the age t th (z, p) in Equation (1), we need an 
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expression for H(z, p). As we have seen, ACDM assumes specific constituents in the 
density, written as p = p r + p m + Pde■ These densities are often written in terms of today’s 
critical density, p c = 3 c 2 H 2 /8ttG, represented as f) m = p m /p c , fl r = p r /p c , and fide = Pde/Pc- 
Hq is the Hubble constant. In a flat universe with zero spatial curvature, the total scaled 
energy density is fi = fi m + fi r + fide = 1- When dark energy is included with an unknown 
equation-of-state, pd e = WdePde, the most general expression for the Hubble parameter is 

H(z, p) = Hq [fi m (l + Z ) 3 + fi r (l + zY + fife(l + Z ) 2 + fide(l + ^) 3 b +,u de)j !■ (2) 

where fi^ is defined similarly to fi m and represents the spatial curvature of the Universe. 
Of course, fi r (~ 5 x 1CU 5 ) is known from the current temperature (~ 2.725 K) of the 
CMB, and the value of Hq. In addition, we assume a flat ACDM cosmology with fi*, = 0, 
for which fide = 1 — fi m — fir, thus avoiding the introduction of fide as an additional free 
parameter. So for the basic ACDM model, p includes three free parameters: fi m , Wd e , and 
H 0 , though to be as favorable as possible to this model, we will also assume that dark 
energy is a cosmological constant, with Wde = —1, leaving only two adjustable parameters. 
However, when we consider the concordance model with the assumption of prior parameter 
values, the fits have zero degrees of freedom from the cosmology itself. 

The Rh = ct Universe is a flat Friedmann-Robertson-Walker (FRW) cosmology 
that strictly adheres to the constraints imposed by the simultaneous application of the 
cosmological principle and Weyl’s postulate (Mclia 2012b; Mclia & Shevchuk 2012). When 
these ingredients are applied to the cosmological expansion, the gravitational horizon 
i?h = c/H must always be equal to ct. This cosmology is therefore very simple, because 
a(t) oc t, which also means that 1 + z = 1/t, with the (standard) normalization that 
a(to) = 1. Therefore, in the Rh = ct Universe, we have the straightforward scaling 

H(z,p) = (l + z)H 0 . (3) 

Notice, in particular, that the expansion rate H(z) in this model has only one free 



parameter, i.e., p is H 0 . From Equation (3), the age of the = ct Universe at redshift z is 


simply 


t th (z,H 0 ) 


1 

(l + z)H 0 ' 


( 4 ) 


To carry out the age-redshift analysis of ACDM and R h = ct , we will first attempt to 
fit the ages of 32 old passive galaxies distributed over the redshift interval 0.117 < z < 1.845 
(Simon et ah 2005), listed in Table 1 of Samushia et al. (2010), assuming a uniform 
value of the time delay r for every galaxy. Following these authors, we will also assume 
a 12% one-standard deviation uncertainty on the age measurements (Dantas et ah 2009, 
2011; Samushia et ah 2010). The total sample is composed of three sub-samples: 10 held 
early-type galaxies from Treu et ah (1999, 2001, 2002), whose ages were obtained by using 
the SPEED models of Nolan et ah (2001) and Jimenez et ah (2004); 20 red galaxies from 
the publicly released Gemini Deep Deep Survey (GDDS), whose integrated light is fully 
dominated by evolved stars (Abraham et ah 2004; McCarthy et ah 2004); and the 2 radio 
galaxies hBDS 53W091 and LBDS 53W069 (Dunlop et ah 1996; Spinrad et ah 1997). 
These data were first collated by Jimenez et ah (2003). 


3. Optimization of the Model Parameters Using a Uniform r 

In subsequent sections of this paper, we will study the impact of a distributed 
incubation time on the overall fits to the data. However, in this first simple approach (see, 
e.g., Dantas et ah 2009, 2011; Samushia et ah 2010) we will assume an average delay factor 
(r) and use it uniformally for every galaxy, so that we may compare our results to those of 
previous work. 

For each model, we optimize the fit by finding the set of parameters (p) that minimize 
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the y 2 , using the statistic 

,2 / v y'[t th (zi,P)r!t G (zi)-(T )] 2 
Xage l T > PJ rr 2 

i=l 

= A — 2 * (t) * B + (t) 2 * C , (5) 

where A = E [t th ( z i, P) - t G {zi )] 2 /cr? Gi , B = E [t th (zi, P) - t G {zi)\ /a? G ., and 
C = E 1 /E c , ■ The dispersions cr tG . represent the uncertainties on the age measurements of 
the sample galaxies. Given the form of Equation (5), we can marginalize (r) by minimizing 
X 2 ag e) which has a minimum at (r) = B/C, with a value x 2 age = A — B 2 /C. Note that 
this procedure allows us to determine the optimized value of (r) along with the best-fit 
parameters of the model being tested. 


3.1. ACDM 

In the concordance ACDM model, the dark-energy equation of state parameter, w^ e , 
is exactly —1. The Universe is flat, fide = 1 — fl m — fir, so there remain only two free 
parameters: fl m and H 0 . Type la SN measurements (see, e.g., Garnavich et al. 1998; 
Perlmutter et al. 1998, 1999; Riess et al. 1998; Schmidt et al. 1998), CMB anisotropy 
data (see, e.g., Ratra et al. 1999; Podariu et al. 2001; Spergel et al. 2003; Komatsu et 
al. 2009, 2011; Hinshaw et al. 2013), and baryon acoustic oscillation (BAO) peak length 
scale estimates (see, e.g., Percival et al. 2007; Gaztahaga et al.2009; Samushia & Ratra 
2009), strongly suggest that we live in a spatially flat, dark energy-dominated universe with 
concordance parameter values fl m ~ 0.27 and Hq a; 70.0 km s -1 Mpc -1 . 

In order to gauge how well ACDM and Rh = ct account for the galaxy age-redshift 
measurements, we will first attempt to fit the data with this concordance model, using 
prior values for all the parameters but one, i.e., the unknown average delay time (r). We 
will improve the fit as much as possible by marginalizing (r), as described above. Fitting 
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the 32 age-redshift measurements with a theoretical t th (z) function using f2 m = 0.27 and 
H 0 = 70.0 km s -1 Mpc -1 , we obtain an optimized delay factor (r) = 1.36 Gyr, and a 
X^ of = 16.17/31 = 0.523, remembering that all of the ACDM parameters are assumed to 
have prior values, except for (r). 



Fig. 1.— l-3cr-constraint contours for the flat ACDM model, using the 32 age-redshift data. 
The cross indicates the best-fit pair (D m , H 0 ) = ( 0.12, 94.3). 

If we relax the priors, and allow both fl m and H 0 to be free parameters, we obtain 
best-fit values (fl m , H 0 ) =(0.12, 94.3 km s -1 Mpc -1 ), as illustrated in Figure 1. The delay 
factor corresponding to this best fit is (r) = 1.62 Gyr, with a y| of = 12.42/29 = 0.428. 
Figure 1 also shows the l-3a constraint contours of the probability function in the Q m -Ho 
plane. Insofar as the ACDM model is concerned, this optimization has improved the 
quality of the fit as gauged by the reduced y^ of , though the parameter values are quite 
different from those of the concordance model. Nonetheless, these two sets of values are still 
marginally consistent with each other because the data are not good enough yet to improve 
the precision with which D m and H 0 are determined. It is also possible that treating (r) as 
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a uniform variable for all galaxies may be over-constraining, but we will relax this condition 
in subsequent sections and consider situations in which r* may be different for each galaxy 
in the sample. We shall see that these optimized values change quantitatively, though the 
qualitative results and conclusions remain the same. The contours in Figure 1 show that at 
the lcr-level, we have 58.5 < H 0 < 127.0 km s^ 1 Mpc^ 1 , and 0.01 < Q m < 0.66. The cross 
indicates the best-fit pair. For the sake of a direct one-on-one comparison between ACDM 
and the Rh = ct Universe, the current status with these data therefore suggests that we 
should use the concordance parameter values, which are supported by many other kinds of 
measurements, as described above. 


3.2. The R h = ct Universe 

Regardless of what constituents may be present in the cosmic fluid, insofar as the 
expansion dynamics is concerned, the R j, = ct Universe always has just one free parameter, 
H 0 . The results of fitting the age-redshift data with this cosmology are shown in Figure 2 
(solid line). We see here that the best fit corresponds to H 0 = 67.2+1'^ km s^ 1 Mpc^ 1 (la). 
The corresponding delay factor is (r) = 2.72 Gyr. With 32 — 2 = 30 degrees of freedom, we 
have x'Lf = 13.05/30 = 0.435. 

To facilitate a direct comparison between ACDM and Rh = ct, we show in Figure 3 
the galaxy ages (i.e., to + (r)), together with the best-fit theoretical curves for the R h = ct 
Universe (with H 0 = 67.2 km s -1 Mpc -1 and (r) = 2.72 Gyr), the concordance model 
(with (r) = 1.36 Gyr and prior values for all the other parameters), and for the optimized 
ACDM model (with H 0 = 94.3 km s -1 Mpc -1 , ff m = 0.12, and (r) = 1.62 Gyr). As 
described above, (r) is the average incubation time or delay factor, which accounts for our 
ignorance concerning the amount of time elapsed since the big bang to the initial formation 
of the object. At the very minimum, (r) must be greater than ~ 300 Myr, this being the 
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time at which Population III stars would have established the necessary conditions for the 
subsequent formation of Population II stars (see, e.g., Melia 2014b and references cited 
therein). The galaxies could not have formed any earlier than this, based on the physics we 
know today. 



Fig. 2.— Constraints on the Hubble constant, Ho, in the context of R h = ct. 

In this figure, we also show to + 300 Myr versus z (circles) to illustrate the minimum 
possible ages the galaxies could have, given what we now know about the formation of 
Population II and Population III stars. We note that the best fit value of (r), in both 
ACDM and R h = ct, is fully consistent with the supposition that all of the galaxies should 
have formed after the transition from Population III to Population II stars at f ~ 300 Myr. 
Strictly based on their x\ Q i values, the concordance ACDM model and the i?h = ct Universe 
appear to fit the passive galaxy age-redshift relationship (i.e., to + (r) versus z) comparably 
well. However, because these models formulate their observables (such as the theoretical 
ages in Equations 1 and 4) differently, and because they do not have the same number of 
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0.0 0.4 0.8 1.2 1.6 2.0 
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Fig. 3.— The complete age-redshift sample ( solid points ), and the best-fit theoretical curves: 
( dot-dashed line ) the concordance model, with its sole optimized parameter (r) = 1.36 Gyr; 
(dashed line ) the standard, flat ACDM cosmology, with optimized parameters fl m = 0.12, 
Hq = 94.3l35;g (lcr) km s -1 Mpc -1 , and (r) = 1.62 Gyr; ( solid line ) the R h = ct Universe, 
with H 0 = 67.2^;o (la) km s -1 Mpc -1 and (r) = 2.72 Gyr. The empty circles show the 
minimum ages the galaxies could have using tc{zi) + 300 Myr. 
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free parameters, a comparison of the likelihoods for either being closer to the ‘true’ model 
must be based on model selection tools, which we discuss in § 5 below. But first, we will 
strengthen this analysis by considering possibly more realistic distributions of the delay 
time r. 


4. Optimization of the Model Parameters Using a Distributed r 

We now relax the constraint that the time delay should have the same value (r) for 
every galaxy, and instead consider two representative distributions: (i) a Gaussian 

P(t) oc exp - ^ , (6) 

where the mean value r c is to be optimized for each theoretical fit, given some dispersion 
ay; and (ii) a top-hat 

{ const (t c — cr T < t < t c -\- ay) 

( 7 ) 

0 (otherwise) , 
with ay now representing the width of the distribution. 

To isolate the various influences as much as possible, we begin with the concordance 
ACDM model, for which Q m = 0.27 and Wd e = —1 (he., dark energy is assumed to be a 
cosmological constant), though we optimize the Hubble constant to maximize the quality of 
the fit. For each distribution P(r) and assumed value of oy, we randomnly assign the time 
delay r* to each galaxy and then find the best-fit values of H 0 and r c by minimizing the y 2 
using the statistic 

2 , \ ^ [f th (*i,p) -t G (Zi) -Ti(T c ,a T )] 2 

Xage'jc , P) = ^ -~2- • ( 8 ) 

i= 1 

The left-hand panels of Figure 4 show the ensuing distributions of r c and H 0 values for 
a Gaussian P(r), and three different assumed dispersions ay. In this case, the optimized 
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H Q (km s ' 1 Mpc 1 ) z 


Fig. 4.— ACDM with a Gaussian distribution of r values. Left-hand panels: fitted values 
of t c and H 0 using the complete age-redshift sample (right-hand panels; solid points). The 
theoretical curves (right-hand panel; dot-dashed curves ) correspond to the parameter values 
that minimize the y 2 (shown on the left). For ACDM, a Gaussian distribution in r results in 
an optimized value of the Hubble constant (~ 77.5 km s -1 Mpc -1 ) only weakly dependent 
on cr r , and a mean delay time r c rs^ 0.86 Gyr. 


age (Gyr) 
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Hubble constant (~ 77.5 km s _1 Mpe -1 ) is effectively independent of oy, while the best-fit 
value of the mean delay time r c is restricted to ~ 0.86 Gyr. Not surprisingly, the scatter 
about the best-fit theoretical curve worsens as oy increases, resulting in larger values of 

y 2 . . 

Assuming a top-hat P(t) with ACDM produces the results shown in Figure 5. The 
best-fit values of Hq and r c are very similar to those associated with Figure 4. In this case, 
both H 0 and r c are effectively independent of the assumed distribution width oy. 

We follow the same procedure for the Rh = ct Universe, first considering a Gaussian 
distribution of r* values (Figure 6), followed by the top-hat distribution (Figure 7). The 
comparison between these two cosmological models may be summarized as follows: the 
best-fit results are very similar for both the Gaussian and top-hat P(t) distributions, for 
the same cosmological model; strictly based on their minimum y 2 values, the concordance 
ACDM model and the i?h = ct Universe appear to fit the passive galaxy age-redshift 
relationship (i.e., to + t* versus z) comparably well, independent of what kind of time-delay 
distribution is assumed. 


5. Model Selection Tools 

Several model selection tools used in cosmology (see, e.g., Melia & Maier 2013, and 
references cited therein) include the Akaike Information Criterion, AIC = y 2 + 2 n, where 
n is the number of free parameters (Liddle 2007), the Kullback Information Criterion, 
KIC = y 2 + 3n (Cavanaugh 2004), and the Bayes Information Criterion, BIC = y 2 + (ln N)n, 
where N is the number of data points (Schwarz 1978). In the case of AIC, with AICq, 


characterizing model A4 a , the unnormalized confidence that this model is true is the Akaike 
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Fig. 5.— Same as Figure 4, except now for the top-hat distribution P(r) given in Equa¬ 
tion (7). In this case, both H 0 and r c are effectively independent of the assumed distribution 
width ay. 
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Fig. 6. 


Same as Figure 4, except now for the R ^ = ct Universe. 
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Best-fit: H =67.0, t =2.72. 
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Fig. 7.— Same as Figure 6, except now for the top-hat distribution P{r) given in Equa¬ 
tion (7). 
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weight exp(—AIC a /2). Model M. a has likelihood 

exp(—AIC q /2) 


P(M a ) = 


exp(—AICi/2) + exp(—AIC 2 /2) 
of being the correct choice in this one-on-one comparison. Thus, the difference 
AAIC = AIC 2 — AICi determines the extent to which A4i is favoured over A4 2 - For 
Kullback and Bayes, the likelihoods are defined analogously. 


( 9 ) 


For the case of the average delay factor (r), with the optimized fits we have reported 
in this paper, our analysis of the age- z shows that the KIC does not favour either = ct 
or the concordance model when we assume prior values for all of its parameters. The 
calculated KIC likelihoods in this case are ~ 51.5% for R\, = ct, versus ~ 48.5% for ACDM. 
However, if we relax some of the priors, and allow both Q m and Ho to be optimized in 
ACDM, then = ct is favoured over the standard model with a likelihood of 66.5% 
versus 33.5% using AIC, « 76.6% versus ~ 23.4% using KIC, and ~ 80.5% versus ~ 19.5% 
using BIC. 


For the distributed time delays discussed in Section 4, the model selection criteria result 
in the likelihoods shown in Table 1. Note that in this case, both the concordance ACDM 
and i?h = ct models have the same free parameters (i.e., Ho and r c ), so the information 
criteria should all provide the same results. For the sake of clarity, we therefore show only 
the AIC results in Table 1, where we see that the AIC does not favour either A h = ct or the 
concordance model, regardless of which distribution is adopted for the incubation time. 


6. Numerical Simulations 

The results of our analysis suggest that the measurement of galaxy ages may be used to 
identify the preferred model in a one-on-one comparison. In using the model selection tools, 
the outcome A = AICi— AIC 2 (and analogously for KIC and BIC) is judged ‘positive’ 
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in the range A = 2 — 6, and ‘strong’ for A > 6. As we have seen, the adoption of prior 
values for the parameters in ACDM produces comparable likelihood outcomes for both 
models, regardless of whether we assume a uniform time delay (r), or a distribution of 
values. Of course, a proper statistical comparison between ACDM and R\, = ct should not 
have to rely on prior values, particularly since the optimized cosmological parameters differ 
from survey to survey. If we don’t assume prior values for the parameters in ACDM, and 
optimize them to produce a best fit to the galaxy-age data, the corresponding A using the 
currently known 32 galaxy ages falls within the ‘positive’ range in favor of R.\ t = ct , though 
not yet the strong one. These results are therefore suggestive, but still not sufficient to 
rule out either model. In this section, we will therefore estimate the sample size required 
to significantly strengthen the evidence in favour of Rh = ct or ACDM, by conservatively 
seeking an outcome even beyond A ~ 11.62, i.e., we will see what is required to produce a 
likelihood ~ 99.7% versus ~ 0.3%, corresponding to 3 a. 

Since the results do not appear to depend strongly on whether one chooses a uniform 
(r) for all the galaxies, or a distribution of individual r t values, we will first use the same 
average value (r) in our simulations, and then discuss how these results would change for 
a distributed incubation time. We will consider two cases: one in which the background 
cosmology is assumed to be ACDM, and a second in which it is R,h = ct, and we will 
attempt to estimate the number of galaxy ages required in each case in order to rule out 
the alternative (incorrect) model at a ~ 99.7% confidence level. The synthetic galaxy ages 
are each characterized by a set of parameters denoted as ( z , t[z]), where t(z) = to + (r). 
We generate the synthetic sample using the following procedure: 

1. Since the current 32 old passively evolving galaxies are distributed over the redshift 
interval 0.117 < z < 1.845, we assign z uniformly between 0.1 and 2.0. 

2. With the mock z, we first infer t(z) from Equations (1) and (4) corresponding either 
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to a flat ACDM cosmology with = 0.27 and Hq = 70 km s~ x Mpc -1 (§ 4.2), or the 
i?h = ct Universe with H 0 = 70 km s -1 Mpc -1 (§ 4.1). We then assign a deviation (At) to 
the t(z) value, i.e., we infer t'(z) from a normal distribution whose center value is t(z) and 
a = 0.35 is its deviation (see Bengaly et al. 2014). The typical value of a = 0.35 is taken 
from the current (observed) sample, which yields a mean and median deviation of a — 0.38 
and 0.33, respectively. 

3. Since the observed error a t is about 12% of the age measurement, we will also assign 
a dispersion cr t = 0.12 t'(z) to the synthetic sample. 

This sequence of steps is repeated for each galaxy in the sample, which is enlarged until 
the likelihood criterion discussed above is reached. As with the real 32-ages sample, we 
optimize the model fits by minimizing the y 2 function y 2 = ^ j [t th (^,p) -t'(zi)] 2 /al^. 
This minimization is equivalent to maximizing the likelihood function £ oc exp (—y 2 /2). 
We employ Markov-chain Monte Carlo techniques. In each Markov chain, we generate 10 5 
samples according to the likelihood function. Then we derive the cosmological parameters 
from a statistical analysis of the sample. 


6.1. Assuming = ct as the Background Cosmology 

We have found that a sample of at least 350 galaxy ages is required in order to rule 
out ACDM at the ~ 99.7% confidence level. The optimized parameters corresponding to 
the best-fit ACDM model for these simulated data are displayed in Figure 8. To allow 
for the greatest flexibility in this fit, we relax the assumption of flatness, and allow f^e 
to be a free parameter, along with Q m . Figure 8 shows the 1-D probability distribution 
for each parameter (D m , Dd e , H 0 ), and 2-D plots of the la and 2 a confidence regions 
for two-parameter combinations. The best-fit values for ACDM using the simulated 
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Table 1. AIC Likelihood Estimation 


Distributed ' 

r a r 

ACDM 

tj 

II 

A 

Gaussian 

0.1 

50% 

50% 


0.3 

47% 

53% 


0.5 

41% 

59% 

Top-hat 

0.3 

51% 

49% 


0.5 

50% 

50% 


0.8 

52% 

48% 



H 0 Q de 


Fig. 8.— The 1-D probability distributions and 2-D regions with the lcr and 2cr contours 
corresponding to the parameters fl m , fide, and H 0 in the best-fit ACDM model, using the 
simulated sample with 350 ages, assuming = ct as the background cosmology. 
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sample with 350 ages in the Rh = ct Universe are f2 m = 0.011, fide = 0.37 +q; 32 (lcr), and 
Hq = 79.9^32v (I 0- ) km s -1 Mpc -1 . Note that the simulated ages provide a good constraint 
on fide; but only a weak one on fl m ; only an upper limit of ~ 0.025 can be set at the lcr 
confidence level. 



Fig. 9.— The 1-D probability distribution for the parameter H 0 in the Rh = ct universe, 
using a sample of 350 ages, simulated with Rh = ct as the background cosmology. The 
assumed value for H 0 in the simulation was H 0 = 70 km s -1 Mpc -1 . 

In Figure 9, we show the corresponding 1-D probability distribution of H 0 for the 
Rh = ct universe. The best-fit value for the simulated sample is H 0 = 70.2+Qcj (lcr) km s -1 
Mpc -1 . The assumed value for Hq in the simulation was Hq = 70 km s -1 Mpc -1 . 

Since the number N of data points in the sample is now much greater than one, 
the most appropriate information criterion to use is the BIC. The logarithmic penalty in 
this model selection tool strongly suppresses overfitting if N is large (the situation we 
have here, which is deep in the asymptotic regime). With N = 350, our analysis of the 
simulated sample shows that the BIC would favour the Rh = ct Universe over ACDM by an 
overwhelming likelihood of 99.7% versus only 0.3% (i.e., the prescribed 3<r confidence limit). 
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6.2. Assuming ACDM as the Background Cosmology 

In this case, we assume that the background cosmology is ACDM, and seek the 
minimum sample size to rule out R h = ct at the 3a conhdence level. We have found that 
a minimum of 45 galaxy ages are required to achieve this goal. To allow for the greatest 
flexibility in the ACDM fit, here too we relax the assumption of flatness, and allow f2 de to 
be a free parameter, along with Q m . In Figure 10, we show the 1-D probability distribution 
for each parameter (O m , f2 de , H 0 ), and 2-D plots of the 1 a and 2a conhdence regions for 
two-parameter combinations. The best-fit values for ACDM using this simulated sample 
with 45 galaxy ages are f2 m = 0.28^g;^ (la), fide = 0.30, and H 0 = OO.ll^g (la) km s' 1 
Mpc -1 . Note that the simulated ages now give a good constraint on fl m , but only a weak 
one on f2 de ; only an upper limit of ~ 0.83 can be set at the la conhdence level. 

The corresponding 1-D probability distribution of Hq for the Rh = ct universe is shown 
in Figure 11. The best-ht value for the simulated sample is Hq = 85.Ol}'® (la) km s -1 
Mpc -1 . This is similar to that in the standard model, but not exactly the same, reaffirming 
the importance of reducing the data separately for each model being tested. With N = 45, 
our analysis of the simulated sample shows that in this case the BIC would favour ACDM 
over R^ = ct by an overwhelming likelihood of 99.7% versus only 0.3% (i.e., the prescribed 
3a conhdence limit). 

These results were obtained assuming a uniform incubation time (r) throughout the 
mock sample. Of course, if the incubation time is distributed, the corresponding uncertainty 
in its distribution function will contribute to the variance of the “observed” age of the 
Universe. Thus, adding the scatter in the uncertain incubation time r c to the simulations 
would change the constructed sample size required to achieve the 3a results discussed above. 
We have therefore carried out additional simulations using a Gaussian P(t) distribution of 
the incubation time, and an assumed dispersion cr T = 0.3. The corresponding uncertainty 
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Fig. 10.— Same as Figure 8, except now with a flat ACDM as the (assumed) background 
cosmology. The simulated model parameters were f2 m = 0.27 and Hq = 70 km s -1 Mpc -1 . 
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in r contributes to the variance of the “observed” age of the Universe. From these results, 
we estimate that a sample of about 55 galaxy ages would be needed to rule out Rh = ct at 
a ~ 99.7% confidence level if the real cosmology were ACDM, while a sample of at least 
500 ages would be needed to similarly rule out ACDM if the background cosmology were 
instead Rh = ct. 


7. Discussion and Conclusions 

In this paper, we have used the sample of high-redshift galaxies with measured 
ages to compare the predictions of several cosmological models. We have individually 
optimized the parameters in each case by minimizing the y 2 statistic. Using a sample of 
32 passively evolving galaxies distributed over the redshift interval 0.117 < z < 1.845, we 
have demonstrated how these age-redshift data can constrain parameters, such as H 0 and 
Q m . For ACDM, these data are not good enough to improve upon the concordance values 
yet, but are approaching the probative levels seen with currently available gamma-ray burst 
luminosity data (Wei et al. 2013), strong gravitational-lensing measurements (see, e.g., 
Suyu et al. 2013), and measurements of the Hubble parameter as a function of redshift 
(Melia & Maier 2013). 

Based solely on these 32 passively evolving galaxies, a comparison of the y 2 iof for the 
Rh = ct Universe and the concordance ACDM model shows that the age-redshift data do 
not yet favour either model. The Rh = ct Universe fits the data with y^ of = 0.435 for a 
Hubble constant H 0 = 67.2tf^ km s -1 Mpc -1 and a average delay time (r) = 2.72 Gyr. By 
comparison, the concordance model fits these same data with a reduced y)j of = 0.523, with 
a delay time (r) = 1.36 Gyr. Both are consistent with the view that none of these galaxies 
should have started forming prior to the transition from Population III to Population II 
stars at ~ 300 Myr. However, if we relax some of the priors, and allow both fl m and Hq 



to be optimized in ACDM, we obtain best-fit values f2 m = 0.121 ojli and H 0 = 94.3^g5;g 
km s _1 Mpc” 1 . The delay factor corresponding to this best fit is (r) = 1.62 Gyr, with a 
Xdof = 0.428. The current sample favours = ct over the standard model with a likelihood 
of « 66.5% - 80.5% versus « 19.5% - 33.5%. 

We also analyzed the age-redshift relationship in cases where the delay factor r may 
be different from galaxy to galaxy, and considered two representative distributions: (i) a 
Gaussian; and (ii) a top-hat. We found that the optimized cosmological parameters change 
quantitatively, though the qualitative results and conclusions remain the same, independent 
of what kind of the distribution one assumes for r. Though one does not in reality expect 
the delay factor to be uniform, the fact that its distribution does not significantly affect 
the results can be useful for a qualitative assessment of the data, ft also suggests that the 
outcome of our analysis is insensitive to the underlying assumptions we have made. 

But though galaxy age estimates currently tend to slightly favour R\, = ct over ACDM, 
the known sample of such measurements is still too small for us to completely rule out 
either model. We have therefore considered two synthetic samples with characteristics 
similar to those of the 32 known age measurements, one based on a ACDM background 
cosmology, the other on 7% = ct. From the analysis of these simulated ages, we have 
estimated that a sample of about 45 — 55 galaxy ages would be needed to rule out 7% = ct 
at a ~ 99.7% confidence level if the real cosmology were in fact ACDM, while a sample of 
350 — 500 ages would be needed to similarly rule out ACDM if the background cosmology 
were instead = ct. These ranges allow for the possible contribution of an uncertainty in 
r to the variance of the observed age of the Universe at each redshift. The difference in 
required sample size is due to ACDM’s greater flexibility in fitting the data, since it has a 
larger number of free parameters. 


Both the Gaussian and Top-hat distributions that we have incorporated into this 
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study have assumed that the mean and scatter of the incubation time are constant with 
redshift. However, it would not be unreasonable to suppose that these quantities could have 
a systematic dependence on z. To examine how the results might change in this case, we 
have therefore also analyzed the real data using a Gaussian distribution 


P(t) oc exp 


T ~ T c • (1 + Zj) a } 2 
2 [oy • (1 + Zi) a ] 2 ’ 


( 10 ) 


with a T = 0.3 and, for simplicity, a — 1. Figure 12(a) shows the corresponding distributions 
of r c and Hq for the concordance ACDM model, with best-fit values ( Hq,t c ) = (82.5,0.26). 
The analogous distributions for the R h = ct Universe are shown in Figure 12(b). In this 
case, the best fit corresponds to (. H 0 ,t c ) = (80.0,0.73). The added redshift dependence has 
not changed the result that both models fit the passive galaxy age-redshift relationship 
comparably well, based solely on their reduced y 2, s. Note that in this case, both the 
concordance ACDM and R,\ t = ct models have the same free parameters (i.e., (H 0 and r c ), 
so the information criteria should all provide the same results. Therefore, we show only the 
AIC results here. We find that the AIC does not favour either = ct or the concordance 
ACDM model, with relatively likelihoods of ss 54% versus ~ 46%. Note, however, that the 
best-fit value of t c in ACDM does not appear to be consistent with the supposition that all 
of the galaxies should have formed after the transition from Population III to Population II 
stars at f ~ 300 Myr. 


An additional limitation of this type of work is the degree of uncertainty in the 
galaxy-age measurement itself. It is difficult to precisely constrain stellar ages for 
systems that are spatially resolved using stellar evolution models. We may be grossly 
underestimating how uncertain the age measurements of distant galaxies are. One ought 
to acknowledge this possibility and consider its impact on cosmological inferences. For 
example, redoing our analysis using a Gaussian P(r) and an assumed dispersion a T = 0.3, 
but now with an additional 24% uncertainty on the age measurements, (i.e., twice as big 
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Fig. 11.— Same as Figure 9, except now with ACDM as the (assumed) background cosmol¬ 
ogy- 
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Fig. 12.— (a): ACDM with a Gaussian distribution for r and a redshift dependent dispersion 
and mean value (see text), (b): Same as (a), except now for the R h = ct Universe. 









31 


as the value quoted in Dantas et al. 2009, 2011, and Samushia et al. 2010), produces the 
results shown in Figure 13. Panel (a) in this plot shows the corresponding distributions of 
r c and H 0 in the concordance ACDM model, with best-fit values ( H 0 , r c ) = (77.5, 0.86). Not 
surprisingly, a comparison of Figure 13(a) with the left-middle panel in Figure 4 shows that, 
as the uncertainty in the age measurement increases, the constraints on model parameters 
weaken; nonetheless, the best-fit values of H 0 and r c are more or less the same. 

The comparison between ACDM (Figure 13a) and R h = ct (Figure 13b) may be 
summarized as follows: the best-fit results are more or less the same for both the 12% and 
24% uncertainties. Based solely on their minimum y 2 values, both models fit the passive 
galaxy age-redshift relationship comparably well. The AIC does not favour either = ct 
or the concordance ACDM model, regardless of how uncertain the galaxy-age measurements 
are, with relative likelihoods of ~ 51% versus ~ 49%. 




Fig. 13.— (a): ACDM with a Gaussian distribution of r and an (assumed) uncertainty of 
24% in the age measurement, (b): Same as (a), except now for the R h = ct Universe. 
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